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A binary system composed of an oscillating and rotating coplanar dusty 
disk and a point mass is considered. The conservative dynamics is treated on 
the Newtonian level. The effects of gravitational radiation reaction and wave 
emission are studied to leading quadrupole order. The related waveforms are 
given. The dynamical evolution of the system is determined semi-analytically 
exploiting the Hamiltonian equations of motion which comprise the effects 
both of the Newtonian tidal interaction and the radiation reaction on the 
motion of the binary system in elliptic orbits. Tidal resonance effects between 
orbital and oscillatory motions are considered in the presence of radiation 
damping. 

1. Introduction 

Compact binaries are among the most promising sources for the gravitational waves 
that should be detected by gravitational wave observatories such as GEO600, LIGO, 
TAMA, and VIRGO on Earth, or LISA in space. On the other hand, it is well known 
that oscillating compact objects also emit gravitational waves that should be detectable 
by gravitational wave observatories. Though the ring-down phases of merged binaries 
or collapsed objects are the strongest sources for gravitational waves from oscillating 
objects, oscillations of the components of binary systems generated by tidal interaction 
during their inspiralling phase should induce measurable characteristic modifications in 
the gravitational waveforms compared with those from binaries with point-like compo- 
nents. Particularly in the leading resonance cases where the orbital frequency is a factor 
of two or three smaller than the oscillation frequencies, the internal motion of the bina- 
ries is expected to be strongly excited. 
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The treatment of the motion of compact binaries within the full Einstein theory of 
gravity is a demanding challenge. Much effort has been invested in this challenge on the 
numerical simulation side because an exact analytic solution cannot be expected to be 
achieveable. However, a fully numerically simulation is still far ahead. 

The aim of the present paper is the development of a quite simple approximate model 
to the full Einstein theory which shows typical features of inspiralling compact binaries 
and which is fully analytic apart from the numerical integration of a system of ordinary 
differential equations. This model of a binary system is given by a point-like object and 
an oscillating, rotating dusty disk, where the orbital plane and the plane of the disk 
coincide. As there is no pressure involved, an extreme simplification occurs for the field 
equations as well as for the equations of motion. The dynamics is treated at the level of 
Newtonian gravity augmented by the gravitational radiation reaction dynamics, which 
occurs in the full Einstein theory to leading approximation at the order l/c" (c denotes 
the speed of light) which is, counted in powers of 2.5 orders beyond the Newtonian 
order of 1 /c^ . Within this model we further assume that the most general motion of the 
dusty disk, apart from its orbital motion, is rotation and oscillation, i.e. we do not allow 
the disk to leave this configuration. This means that the tidal force is considered as a 
perturbation only. The same is true for the radiation reaction force which is of 2.5pN 
(post-Newtonian) order. 

The Newtonian oscillating rotating disk of dust is interesting for it's own reasons. 
While stellar oscillations are usually treated as small perturbations, the oscillations of 
the disk can be arbitrarily large, they are bulk oscillations. Moreover, the equation of 
motion for the disk, which can be derived from the equations of hydrodynamics, exhibits 
a surprising similarity to the well-known Keplerian problem. Exploiting this analogy, it 
is possible to write down a parametrized solution for the disk's oscillations jlj. 

Having fully specified our model, let us shortly summarize the present knowledge of 
the analytical motion of binary systems in the Einstein theory of gravity now. Regarding 
point-like binaries without proper spins the conservative motion is known fully analyti- 
cally up to the 3pN order, where the 2.5pN order is neglected j2I; the motion of a single 
oscillating and rotating dusty disk in it's rest frame is known fully analytically to IpN 
order [Hj. If the spins are included into the binary motion, in particular, if the spin ori- 
entations are parallel or anti-parallel to the orbital angular momentum as in our model, 
the analytical solution is known up to 2pN order |^. Taking into account Newtonian 
tidal interaction in close binary systems, several semi-analytical investigations to the 
evolution of binaries have been undertaken [2] , [S] , [Zj • In |H! and jHj the semi-analytical 
solutions of close Newtonian binaries systems have been obtained, including 2.5pN grav- 
itational damping effects. 

The paper is organized as follows: In section II we derive the Hamiltonian formalism for 
the Newtonian disk, starting from the equations of hydrodynamics. We will not consider 
quadrupole radiation damping as a secular effect, but include radiation reaction terms in 
the Hamiltonian equations of motion. To that end we derive the leading order radiation 
reaction Hamiltonian. Then we calculate the leading order gravitional waveform for the 
isolated Newtonian disk. 

In section III we give a short review of orbital Hamiltonian formalism and calculate 
the contributions to the leading order gravitational waveform of the orbital motion. 
In section IV we focus on the binary system. The tidal interaction Hamiltonian is 
calculated explicitly and given in form of a series expansion. In a next step we derive 



the radiation reaction Hamiltonian for the coupled binary system. The leading order 
Hamiltonian equations of motion, governing the dynamics of the system, are solved 
numerically for different initial conditions. Since the contribution of the interaction 
of the disk quadrupole with the orbital quadrupole is much smaller than all the other 
contributions it can be neglected in most applications. On the other hand, the influence 
of the tidal interaction potential on the dynamics of the binary is not negligible because 
the exchange of orbital energy with oscillation energy can be quite strongly. 

Taking into account leading order gravitational damping, orbital energy is dragged 
from the orbit via tidal interaction. In other words, tidal interaction speeds up the 
inspiral process of our binary system. 

The energy transfer is strongest in tidal resonance, where the oscillation frequency is 
n times the orbital frequency (n ist a small integer number) . We study the effect of tidal 
interaction on the orbit and on the leading order gravitational waveform. In section V 
we shall discuss our results and give an outlook on future projects. 

2. Oscillations of a rotating disk of dust - Hamiltonian 
fomulation 

One of the simplest axisymmetric systems that exhibit the essential new features of gen- 
eral relativity is a thin disk of dust. During their evolution these disks emit gravitational 
radiation and may eventually form a black hole. No wonder that general relativistic disks 
of dust have been subjects of interest for many years. Bardeen and Wagoner [Oj where 
the first to study relativistic stationary rotating disks. They succeeded to find a solution 
for a rigidly rotating disk in form of a series expansion. Unfortunately this solution 
can be obtained only numerically beyond 2pN approximation. In terms of ultraelliptic 
functions, Neugebauer and Meinel JOl found an exact solution for a rigidly rotating disk 
in full general relativity using inverse scattering methods. However, these solutions are 
restricted to rigidly rotating stationary disks. Nonstationary, rigidly rotating disks of 
dust were first studied by Hunter p. These disks are particularly interesting, as they 
are sources of gravitational waves. In the present paper we will focus on oscillating 
MacLaurin disks. These are infinitely thin, rotating disks of dust which are obtained as 
a limiting case of MacLaurin spheroids. For these configurations Kley and Schafer |H] 
succeeded to derive a fully analytically solution up to first post-Newtonian approxima- 
tion in a parametrized form. That solution exhibits a surprising similarity to the IpN 
parametrized solution of the Keplerian problem. 

However, the study of a weakly relativistic disk as a component of a binary system is 
beyond the scope of this paper. We shall restrict ourselves on leading order results, thus 
considering only the Newtonian disk. 

Let us start by briefly reviewing the known analytic solution for the oscillating disk, 
obtained from the equation of hydrodynamics for a pressureless system. Since we are 
interested in the leading order gravitational waveform emitted by the oscillating rotat- 
ing disk we shall also discuss the related 2.5pN radiation reaction. Note that we do not 
treat quadrupole radiation damping as a secular effect. The corresponding radiation 
reaction terms will be included into the Hamiltonian equations of motion. We shall 
solve the Hamiltonian equations of motion for the isolated Newtonian disk numerically 
and calculate the leading order gravitational waveform in terms of canonical conjugate 
variables. 



2.1. The Newtonian disk revisited 



The Newtonian dynamics of MacLaurin disks has been treated by several authors pQ, 
jllj , |3j , starting with the equations of hydrodynamics in Eulerian form. Let us briefly 
review this standard approach before we consider the disk in Hamiltonian formalism. 
For an axisymmetric, pressureless configuration confined to the z = plane the equations 
of hydrodynamics and mass conservation in cylindrical coordinates (r, ip) read 

as as S 9 , , 

ot or r or 
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Here S(r, t) denotes the surface density and Vr and are radial and azimuthal velocities, 
respectively. The gravitational potential U is the solution to Poisson's equation 

A;7 = -4^GS(r,t)(5(z). (2) 

For S, fr, % shall choose an ansatz where the surface density is that of a MacLaurin 
disk, 



T.{r,t)=a{t)^l- Vrir, t) = -rfit), v^ir,t) = r^it), (3) 

where cr(t) is the surface density at the center and r^(t) denotes the disk radius at a 
given time. The disk is rigidly rotating with angular velocity Using this ansatz, it 

is possible to solve analytically. To that end we note that the Poisson equation (|2|) 
can be solved explicitly inside the disk. Inserting Q we find 



U{r,t)=Ga{t)rd{t)^ 



rdit) 



ir<rd,z = 0). (4) 



Substituting equations © and Q into (0) we can derive a system of four differential 
equations for the unknown variables r^, cr, / and ^l: 

a - 2fa =0, 
rd + frd =0, 

n - 2fn =0, (5) 
-f + f + ^ =0. 

This system is integrated with respect to the initial conditions /(O) = and /'^(O) = Rd- 
Combining the first two equations, one finds 

^(tydit) = ^Md, (6) 

where is the total Newtonian mass of the disk. The third equation in Q gives 

nityjit) = QoRl no = nio) (?) 



which is constant at the Newtonian level. The remaining equation containing r^i reads 

2C 



where 

C := Ga{t)Ta{tf- = —GM^ (9) 

and 

:= 2CeRd = ^iRi (10) 

is defined by as follows, 

= (11) 

Equation (jSJ exhibits a surprising similarity to the radial equation in the Keplerian prob- 
lem. Exploiting this similarity, the solution to equation Q can be given in parametrized 
form, 

27r 

f d = o,d{)- — e cos (f), —t = u — esmu — TT, (12) 
where the definitions 

Rd 



dd 



l + e 



and e = 1 — hold. For not being confused with the orbital eccentricity e, we shall 
speak of e as the ellipticity of the disk's motion. 



2.2. Hamiltonian formulation - Newtonian level 

After this short review we shall now derive the Hamiltonian formulation of the disk's 
dynamics. On the Newtonian level this can be easily achieved by exploiting the analogy 
of (jEJ to the radial equation of the Keplerian problem. The Lagrangian, leading to Eq. 
©, reads 



L 



a 



2 rd 



a, (3 and 7 being coefficients which have to be determined. The Euler-Lagrangian equa- 
tion corresponding to this ansatz reads 



d dL 
dt dr'd 



dL_ 
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Comparing this with Eq. (jSJ, we have to identify 



(3 = 2C, 7 = -y = -CRdf- 



Substituting the canonical conjugate momentum to r^^, 




dL 

the Newtonian Hamiltonian reads 

^r^ = ^-«(f +ik) = -b.^+^i-^> (13) 

where we have identified = 2a^CRd^^- To fix the value of a we require that H^^^^ 
should be the conserved energy of the disk. The later is given by 

In particular, this must hold at t = 0, where by virtue of the initial condition r^(0) = Rd 
and rd(0) = 0, the momentum conjugate has to fullfil Pr(0) = 0. Inserting this into Eq. 
([T3|l. we end up with 

, a = -Md. 14 

2.3. Radiation reaction part 

It is well known that the leading order part of the gravitational radiation is of the 
quadrupolar type. In the case of continuous matter the related radiation reaction Hamil- 
tonian is given by jl2j 

2C f TT-TT- 1 

Hreacit) = ^Qiji*) J d?x[^ + ^diU.djU,), (15) 

where 

Qij := J <fxp{x'x^ - (16) 

is the Newtonian mass-quadrupole tensor of the system. At Newtonian order we may 
identify = p, vr, = v^p, and = U. Using the continuity equation for the mass and 
the equations of motion for the matter, the following equation, relating an integral of 
non-compact support to an integral of compact support, can be shown to hold. 
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d'xipvV + -^diUdjU) = ^Q,j. (17) 



Thus we are allowed to write the reaction Hamiltonian of our disk as 

H^tl{t) = ^Q[f{t)Q[f{p,q), (18) 

where Q^j^ {p, q) denotes that Q^'j^ has to be treated as a function of position and mo- 
mentum variables, i.e. making use of the equations of motion the second time deriva- 
tives have to be eliminated. The mass quadrupole tensor of the disk is diagonal, 
Qij^ ~ {Qii^Qii^~'^Qii)- This simplifies the Eq. ((TH|l enormeously. In the appendix 



the disk's quadrupole tensor and it's time derivatives are given explicitly. In particular, 
we get 



1 
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Inserting this into Eq. (|18|) we finally arrive at 
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2.4. Equations of motion 

The total Newtonian plus leading order radiation reaction (non-autonomous) Hamilto- 
nian of the oscillating rotating disk reads 



H. 
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where H^^^ and HrVac are given by (|14() and H2U|) . respectively. The Hamiltonian equa- 
tions of motion are obtained in the usual way as 
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However, some care is needed when calculating the reaction part of the Hamiltonian 
equations of motion. It is obtained by taking the derivatives of Hrtlc with respect to the 
canonical conjugate variables, but one must take 'Q['j\t) as a function of t. Only after 
the differentiation we shall express the third time derivative of the quadrupole tensor 
as a function of Pr and r^. Inserting Eq. H19|) the reaction part of the Hamiltonian 
equations of motion reads 
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The Hamiltonian equations for the oscillating disk, including leading order gravitational 
radiation reaction, read 
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Figure 1: In leading order approximation ij^,^ is the only nonvanishing contribution to 
the gravitational radiation field of the rotating, oscillating disk. In the figure 
idisk is plotted for = 0.9,i2d(0) = 10 (in units of Ga/c^). 



Note that we do not treat gravitational damping as a secular effect. The leading order 
radiation reaction terms enter directly into the equations of motion. It is remarkable 
that, although the total energy of the disk is not conserved, the angular momentum is 
a conserved quantity at 2.5pN as indicated by Eq. (|25|) . This is due to the rotational 
symmetry of the disk and well expected. For a binary point-mass system, however, 
2.5pN radiation reaction destroys orbital momentum conservation, as we shall see in the 
next section. 



2.5. Gravitational waves from oscillating rotating disks 

The 2.5pN radiation reaction terms in Eqs. (|22j) - (|25|) are the leading order dissipative 
terms. Energy is dragged from the oscillating disk by emission of gravitational waves. 
Asymptotically, the gravitational radiation field can be represented by [l,*- 
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where D denotes the source-observer distance, the indices i,j refer to Cartesian coor- 
dinates in the asymptotic space, j^^^'^™ and T^'^'''"^ are the pure-spin tensor-spherical 
harmonics of magnetic and electric type and j'™" and S^"^ are the spherical radiative 
mass and current multipole moments, respectively. The upper pre-index (l) denotes the 
number of time derivatives. In leading order approximation only I = 2 terms contribute 
to the radiation field, i.e. 
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This is the famous quadrupole radiation. 

In many cases it is more convenient to work with STF mass multipole moments Qai- 



They are related to the sphericaUy radiative mass multipole according to, e.g. |14j 
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where 
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The brackets denote the symmetric part. In particular, we get 
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Inserting these expressions, we readily find 



r20 



IGtt 



IBtt / 1 'i /27r 
/22 ^ ^_^_(g^^ _ - 2iQu) = '^\j^{Qll - Q22 - 2iQi2). 



(30) 
(31) 
(32) 

(33) 
(34) 
(35) 



The quadrupole tensor of the disk is diagonal and, moreover, Q^^^ = Q^22 • Hence only 
■^disk contributes to the leading order gravitational radiation field of the disk. In terms 
of rd,Pr,y^,P(p the second time derivative of (^33^ reads 



2Ca 
rd 



Inserting this into Eq. H33() . the leading order time dependence of the gravitational 
waveform of the oscillating rotating disk of dust is given by 



r20 

disk 
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(36) 



In figure ^ the 1^^ component of the leading order gravitational radiation field is given 
for a particular example. 

Exploiting the analogy of the analytic solution Eq. H12j) we may express I'jfgf^ in terms 
of u. To that end we start with Eq. ()19|). where Qi? is given as a function of the 
Newtonian energy and r^. Inserting 



Ed = - — {l + e), 



rrf = arf(l - ecosw) 



Rd 
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into H19() we readily obtain 
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and hence 
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This expression is quite similar to the result given for l'^^^ (see Eq. 



3. The orbital motion reviewed 
3.1. Hamiltonian formalism 

The dynamics of a binary point-mass system, including leading order gravitational ra- 
diation reaction, is well known, and we shall only briefly review the basic points. The 
orbital Hamiltonian reads 



Horb - H^^jj + H^°l^{t). 



(37) 



Calculating the leading order reaction Hamiltonian according to Eq. (|15jl we arrive 
at [13 
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and the Newtonian Hamiltonian is given by 
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where A4 is the total and /i the reduced mass, respectively. The time evolution of the 
system is governed by the Hamiltonian equations of motion. 
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Note that Q^fjit) has to be taken as a function of time. Only after the differentia- 
tion it shall be expressed as a function of PR,P$,i? and In polar coordinates, the 
Hamiltonian equations of the orbital motion read 
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where v = fj,/M. 

For numerical calculations we shall apply rescaled variables. Defining X,P'j^,P^ and r 

by 
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we arrive at 
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3.2. Leading order gravitational waveforms 

The leading order gravitational radiation field is given in Eq. (|27j) . If we identify the 
orbital plane with the = 7r/2 - plane, Q^^^ and Qga^ vanish and thus only /^o 
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are nonzero. The time derivatives of Q -- as functions of the canonical variables are 
given in the Appendix B. Inserting this in (|35|) and respectively, we find 
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and 
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3.3. Gravitational waveforms - analytical results 

The contributions to the gravitational waveform in Eqs. (|47|) and (|5flj) are given in terms 
of Pjj,P^,<I> and X . These expressions are applied when solving the dynamics fully 
numerically. It is instructive, however, to compare those results with semi-analytical ones 
where the conservative dynamics is solved fully analytically and where the gravitational 
radiation damping is treated as a secular effect. If the two components are well separated, 
the curves turn out to be in good agreement. At distances of around 10 GM./cP', /22 
calculated by the secular approach is in a good agreement to the one calculated by full 



numerical simulation over a few periods only. 

The parametrized solution to the Newtonian problem reads 



X = a'(l — e cos n), 
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-t = n- t = u — e cos u. 



Here u is the mean anomaly, the period P is given by 
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and the azimuthal angle <I> is related to u by 
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From Eq. H52() we readily obtain 

6 



G7W 



1 + e ^ n 

tan 

I ■ 2 



(51) 



(52) 



n\/l — 



n = r 

1 — e cos u (1 — e cos 



(53) 



The Newtonian energy of the orbit \s E = fic^E' with E' = Using this relation, it 

is now easy to derive from the eqs. and (1^5]) the leading order contributions jl4j 
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where A{u) = 1 — ecosn. In particular, I^^ vanishes identically for circular orbits. The 
luminosity of the gravitational wave emission in leading order approximation is given by 
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Averaging C over the orbital period gives the leading order energy loss 
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while the decay of the orbital eccentricity is given by 
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4. The binary system 

In previous chapters we have considered the Newtonian oscillating rotating disk and the 
orbital motion separately. In the following we shall focus on a binary system composed 
of a point-like object and a dusty disk. Due to tidal interaction neither the internal 
dynamics of the disk nor the orbital motion are completely independent. In fact, the 
influence of tidal coupling on the form of the orbit and on the inspiral process can be 
rather strong, as we will see later on. For simplicity we assume the orbital plane to be 
the z = plane, i.e. the plane of the disk. 



4.1. Tidal coupling on the Newtonian level 

The Newtonian interaction potential for a general matter distribution reads 

G ff p{v)p{i) 3 3 
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Not taking into account self- interaction terms, we are left with 

= U!i^ + C/iS" = -g/ f!f£f Mr. (57) 

One of the two volume integrals is trivial since the point mass density is given by Pp{r) = 
Mp6{i — R). The remaining integral reads 

Uint = -GMpl ^^dV (58) 

Inserting the density of the disk p^ = 5{z)T,{r,t), where the surface density S is defined 
by Eq. (jSJ and where the origin of the coordinate system has been shifted to the center 
of the disk, we end up with 
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where J is defined as 
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Unfortunately, this integral seems not to be solvable analytically for R > (in case of 
R < Td, see e.g. Kley and Schafer [3]). Expanding the integrand in terms of Legendre 
polynoms, we obtain 
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Defining 
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we end up with 



00 ^ /'27r 



The integral J^^ Pi{cos(f)d(f is only nontrivial for even numbers of /. With the substi- 
tutions y = cos ip, dip = the case I = 2n can be rewritten as 



P2n{y) 



/ P2n(cOS99')V = 2 / 
Jo J-1 \J\ 



This is a special case of the more general relation 
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which holds for 5RA > Taking A = ^, /U = 0, u = 2n we arrive at 
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Inserting this into Eq. (jHTjl the tidal potential Eq. (|59|l reads 
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Wi := I even, (63) 
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and Wi = for odd values of /. It is easy to see that the / = term in Eq. ()62p gives 
just the point-mass interaction, 

Uo = "i R r(i)2r(5) = = ('"^ 

where we used r(^) = and r(|) = j\/tt. The tidal interaction potential can be 
immediately read off from Eq. (|62l) . 

1=2 

As expected, the first nonvanishing contribution comes from the quadrupole term jl7j . 
Correspondingly, we expect the tidal interaction to introduce a periastron shift into the 
Newtonian orbit (see e.g. fTF). 



4.2. Radiation reaction Hamiltonian and Hamiltonian equations of motion 



Though it is tempting to take the leading order radiation reaction Hamiltonian of the 
binary system as the superposition of the radiation reaction Hamiltonians of the isolated 
disk and pure orbital motion, respectively, a careful analysis shows that coupling between 
orbital and disk quadrupole moment already enters at leading order. This can be easily 
seen from the general 2.5 pN reaction Hamiltonian of a fluid configuration 
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(66) 



Here qIj^ ^-iid Q-°^ are the STF mass quadrupole moments of the disk and the orbit, 
respectively, R^,Pi denote the canonical conjugate variables of the orbital motion and 
x'^,Pi are the canonical conjugate variables of the disk relative to the disk's center. With 
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of the binary system is 



HreLc and Hreac given by the eqs. ([20]) and (|38|) the leading order reaction Hamiltonian 
Hreac{t)=H(iUt)+H(ll{t) + ^[Q[f{t)Q^^^ . (67) 



The last two terms describe the coupling between disk and orbital quadrupole moments. 
Finally, using the diagonal structure of the disk's mass quadrupole tensor we arrive at 
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(68) 



As before, the radiation reaction parts of the Hamiltonian equations of motion are ob- 
tained by differentiating Hreac with respect to the generalized coordinates and momenta. 
Thus the reaction part of the equations of motion reads. 
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Note that the couphng of disk and orbital mass quadrupole moments is aheady present 
at leading order. However, numerically these terms turn out to be very small even in 
close orbits and thus can be neglected in most scenarios. 



4.3. Hamiltonian equations of motion 

The total Hamiltonian of the binary takes the form 

Htot 
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where 
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and where Hreac is given by Eq. l|H5|). The radiation reactive part of the equations of 
motion has been calculated in the previous section and is given by the Eqs. H69() - H76|) . 
The dynamics of the binary, including tidal coupling and leading order gravitational 
radiation reaction effects, is described by the following set of first order differential 
equations: 
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In order to solve these equations numerically we shall introduce scaled variables. We 
choose the same scaling as in Eq. (|42() for the disk variables 
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Substituting this into the Eqs. 1)811 - 1881) and defining 

C 



A :-- 



_ 3iTMd 
GM ~ T AT' 



B 



a _2 M 
/i ~ 5 iWp 



we finally arrive at 
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4.3.1. Initial conditions 

Both, the dynamics of the binary as well as the leading order gravitational waveform 
depend on the choice of initial conditions. Various initial conditions correspond to 
various initial phase differences between orbital and disk angular variables. Let t = be 
the time of the periastron passage where by convention $(0) = 0. The initial values for 
the orbital variables shall read 



X(0) =a'(l-e), $(0) = 0, P4(0) = V«'(l " e^), P'r 
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(98) 



To fix the initial conditions for the disk variables, we require 
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where e is the elhpticity of the disk at i = and is the maximal disk radius. Note 
that if Xrf(O) = Xfi the momentum conjugate p^(0) vanishes. The momentum conjugate 
to 9J is a constant as can be seen from the equations of motion. It is given by 



pI = 2a'CR,e = ^^^^P'^ =^ P'^m = By/2AX,{l-e). 
The initial value of p'^ for a given Xd{0) can be easily derived from Eq. ()99() . 

p;2(o) = 2B^A 

5. Discussion 

In the last decade several authors investigated tidal and tidal-resonant effects in binary 
systems. Close binary neutron stars were studied by e.g., Schafer and Kokkotas and 
Lai and Ho jSi • Recently, resonant excitations of white dwarf oscillations in a compact 
binary were investigated by Rathore et al. ^H]- In these papers the stellar oscillations 
were treated numerically. In our model there is only one oscillation mode, but it is 
known analytically. Although our model is much simpler than the binary neutron stars 
of jn])|Hl or the white dwarf-compact object binary of it exhibits all essential rela- 
tivistic features of a compact binary system to leading order approximation. Moreover, 
we were able to give the Hamiltonian equations of motion in analytic form. 
In their recent paper Rathore et. al. investigated a white dwarf-compact object system 
with non-rotating white dwarf. Gravitational damping was taken into account as a sec- 
ular effect and the change of the orbital motion due to tidal interaction was discussed 
only qualitatively. The symmetry of our disk-compact object binary enabled us to study 
the effects, in particular the influence of tidal interaction on the actual orbit, in great 
detail (see Fig. Note that while the white dwarf of Rathore et al. was nonrotating, 
rotation is essential in our case to stabilize the dusty disk. The single eigenmode of the 
oscillating disk corresponds to the f-mode of Rathore et al. 

In our paper we studied the oscillations of an isolated rotating Newtonian disk of dust in 
Hamiltonian formalism. Thin oscillating disks of dust with a diameter of a few hundred 
kilometers could eventually form after the merging of two white dwarfs. The similarity 
of the disk's equation of motion to the well-known Keplerian problem is remarkable. 
Exploiting this analogy one can easily write down the analytic solution for the disk's 
oscillations. This allows to study the effects of tidal interaction and it's influence on 
the actual gravitational waveform in great detail (see Figs. 131415161715)) . For a suitable 
choice of parameters (i.e. the initial major semi-axis not being too small) the evolution 
of the system can be followed over hundreds or even thousands of orbital periods with 
considerable accuracy (see e.g. Fig. EllHl)- Our model might thus be a good test to run 
data analysis codes that intend to extract physical parameters of the binary system from 
a given gravitational waveform. 

A particular feature of the tidally coupled binary system is the so-called tidal resonance 
(see Fig. IS)). In Newtonian theory the tidal coupling between the internal dynamics of 
a star and the orbital dynamics becomes particularly strong if one of the star's eigen- 
frequencies is n times the orbital frequency. Tidal interaction and in particular tidal 
resonances may speed up the inspiral process by a considerable amount (see Figs. ISJand 
121). In fact, for a realistic choice of parameters (i.e. of and Mp) the 2:1 resonance 



l + e 
' Xd 



+ 



p'^m 



+ 



B^2AXd{l-e), 
2 

Xd,0 



(100) 
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leads to an immediate destruction of the disk-compact object binary. At this point, the 
analysis breaks down. Higher order resonances, starting with the 3:1 one, can be treated 
without problems for a suitable choice of parameters. 

Some remarks should be made on the validity of our model. It relies on the assumption 
that the disk is not getting disrupted and that the change in the disk's density does 
result from the oscillation of the disk only. In particular we do not take into account 
mass accretion on the compact companion. Thus, even for tidal resonance, the maximal 
radius of the disk has to be located far inside it's Roche lobe. In other words, the max- 
imal disk radius must be considerably smaller than the periastron distance between the 
two objects. 

In a forthcoming paper we will extend our analysis to a system of two oscillating disks. 
In this case there will occur quadrupole-quadrupole coupling between the disk's mass 
quadrupole moments. These terms are of course small but they might become im- 
portant when extending the analysis to include first order post-Newtonian corrections. 
By investigating post-Newtonian corrections we hope to get a better intuition for a 
semi-analytical treatment of the much more complicated analysis of inspiraling NS-NS 
binaries. 
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A. The quadrupole tensor of the disk and it's derivatives 

To calculate the quadrupole tensor of the disk, we start with the general definition 



Since pd, the density of the disk, is confined to the z = plane and depends on r, the 



radial coordinate, only, the nondiagonal elements of Q^f^ vanish. It is now easy to see 



that 
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Due to the symmetry it is clear that the relations 
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hold. The second and third time derivative of can be calculated using the equation 
of motion Eq. With 
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B. Time derivatives of the orbital quadrupole tensor 

By definition, the z = plane is the orbital plane. The Qfg** and Q2s^ components of 
the orbital quadrupole tensor thus vanish and we are left with 



Qff = ^[1 + 3008(2$)], 



Q22 

Q33 



6 



6 



[l-3cos(2$)], 



W12 



sin(2$). 



3 ' 2 

The time derivatives are calculated using the (Newtonian) equations of motion, 
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Figure 2: Tidal interaction and gravitational radiation effects are shown for a slightly 
elliptic (e(0) = 3/10) binary with a'(0) = 19,Xd(0) = Xd = Q,Md = Mp = 
1/2, e = 1/10. At t = the disk to orbital frequency relation is uj^j'^o = 7.05. 
The Newtonian point-particle orbit is given in the upper left figure. Including 
tidal interaction leads to a periastron shift (upper right). The figures on the 
bottom are calculated including leading order gravitational radiation effects in 
the equations of motion. Tidal interaction speeds up the inspiralling process 
as can be seen by comparing the orbit including (right) and excluding (left) 
tidal interaction. 
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Fi gure 3: The influence of the tidal interaction on the IfQi^i component of the binary's 
gravitational radiation field is shown. Real and imaginary parts of Ifg^^l 
given without (left) and including (right) tidal interaction. In order to make 
the influence of tidal interaction on the gravitational waveform more visible, 
the orbital separation is taken to be small (Md = Mp = 1/2, e(0) = 0,a'(0) = 
10, Xd{0) = Xd = 3, e = 1/10). The time is measured in units of Tq, the initial 
orbital period. 




Fi gure 4: I^^ (solid line) and /^^^ (dotted line) for the tidally coupled system (left), 
and the If^^^^i component of the binary's leading order gravitational waveform 
(right). The parameters are the same as in Fig. 3. Although the unperturbed 
orbit is assumed to be circular, /^^^ is not zero due to tidal interaction. 
is dominated by the disk's contribution. The time is measured in units of Tq, 
the initial orbital period. 
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Figure 5: The l'^^ component of the leading order gravitational radiation field of the 
binary for a'(0) = 15, all other parameters are the same as in figures (jSJ, ®. 
On the left, the contributions of the disk is given, the right figure shows Ii^i^i' 
which is dominated by the disk's contribution. The time is measured in units 
of the inital orbital period Tq. 
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Figure 6: Real (solid line) and imaginary (dotted line) part of I^^^g^i for an elliptic orbit 



with e(0) = 3/10, a'(0) = 25,Xd(0) = Xa = 6,e = 1/10, 
Time is measured in fraction of the initial orbital period Tq. 
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Figure 7: /^^ for an elliptic orbit (e(0) = 0.3). In the left graph orbital (dotted line) and 

disk (solid line) contributions arc given separately, the right graph shows the 
■^total P^^^ °f binary's leading order gravitational waveform. The time is 
given in fraction of the initial orbital period Tq. 
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Figure 8: The /^"^ contribution to the leading order gravitational wave radiation is shown 
for a binary with a(0) = 100, e(0) = 0.3, e(0) = 0.1,Xd(0) = Xd = 40, Mp = 
= 1/2. The system starts at cu^/cco = 4.95. The upper left figure shows 
-^Jisfe' '^hils is given in the upper right figure. The total contribution 
is shown in the last figure. The time is measured in units of the initial orbital 



